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We consider a type of intermittent behavior that occurs as the result of the interplay between 
dynamical mechanisms giving rise to type-I intermittency and random dynamics. We analytically 
deduce the laws for the distribution of the laminar phases, with the law for the mean length of the 
laminar phases versus the critical parameter deduced earlier [PRE 62 (2000) 6304] being the corollary 
fact of the developed theory. We find a very good agreement between the theoretical predictions 
and the data obtained by means of both the experimental study and numerical calculations. We 
discuss also how this mechanism is expected to take place in other relevant physical circumstances. 
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Introduction 

Intermittency is known to be an ubiquitous phe- 
nomenon in nonlinear science. Its arousal and main sta- 
tistical properties have been studied and characterized 
already since long time ago, and different types of inter- 
mittency have been classified as types I— III p], B, on-off 
intermittency 0, 0, M , eyelet intermittency [3, d, Q 
and ring intermittency (lQj . 

From the other side, increasing interest has been put 
recently in the study of the constructive role of noise 
and fluctuations in nonlinear dynamical systems. In par- 
ticular, it was discovered that random fluctuations can 
actually induce some deg ree of order in a large variety 
of nonlinear systems ["LlT H2L Qjl, and such phenomena 
were widely observed in relevant physical, chemical and 
biological circumstances [H, [3, 0, [EH • 

There are no doubts that different types of intermittent 
behavior may take place in the presence of noise and fluc- 
tuations in a wide spectrum of systems, including cases 
of practical interest for applications in radio engineering, 
medical, physiological and other applied sciences. It is 
plausible that such an interaction would originate new 
types of dynamics. Therefore, the intermittent behav- 
ior in the presence of noise has been studied by means 
of Fokker-Plank equation [ItJ and adopting renormaliza- 
tion group analysis [Hj], but the characteristic relations 
were obtained only in the subcritical region, where the 
intermittent behavior is observed both in the presence 
of noise and without noise. Recently [l9j], the theoretical 
consideration of the intermittent behavior in the presence 
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of noise has been considered in the supercritical region 
(where intermittency is absent in the absence of noise), 
with the analytical form of the dependence of the mean 
length of the laminar phases versus the critical parameter 
being deduced under the assumption of the fixed reinjec- 
tion probability taken in the form of a (^-function. More- 
over, the found analytical law has been verified by means 
of the experimental observation of the characteristic rela- 
tions of intermittency in the presence of noise [2O, Hl| . At 
the same time, the other important statistical character- 
istic of the intermittent behavior, namely, the distribu- 
tion of the laminar phase lengths, has not been obtained 
hitherto for the supercritical parameter region. 

In this paper we report for the first time the form 
of the distribution of the lengths of the laminar phases 
deduced analytically for the type-I intermittency in the 
presence of noise for the region of the supercritical pa- 
rameter values. The already known dependence of the 
mean length of the laminar phases on the criticality pa- 
rameter |l9|, [20[ follows as a corollary of the carried out 
research. Moreover, we prove that this dependence ob- 
tained in [l9[ under the assumption of the fixed reinjec- 
tion probability taken in the form of delta-function does 
not depend practically on the relaminarization proper- 
ties, and, correspondingly, the obtained expression of the 
mean length of the laminar phases on the criticality pa- 
rameter remains correct for different forms of the reinjec- 
tion probability. The obtained analytical distribution of 
the laminar phase length is verified by means of both nu- 
merical calculations of the model system dynamics and 
experimental observations. 

The structure of the paper is the following. In Sec. [J we 
describe the theory of the type-I intermittency with noise 
and give the theoretical predictions for the distributions 
of the laminar phase length. In Sec. [II] we describe the 
dynamical systems used to illustrate our conclusions. We 
show numerically that our theoretical predictions are ob- 
served in the different nonlinear systems, including cou- 
pled chaotic oscillators near the boundary of the phase 
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synchronization in the case of small detuning of the nat- 
ural frequencies. Finally, in Sec. IIIII we give the descrip- 
tion of the experimental setup for the measurement of the 
characteristics of the type-I intermittency in the presence 
of noise. The final conclusions are given in Sec. IIV1 

I. THE THEORY OF THE TYPE-I 
INTERMITTENCY IN THE PRESENCE OF 
NOISE 

The standard model that is used to study the type-I 
intermittency [l| is the one-parameter quadratic map 

%n+l f(x n ) =X n +X 2 n +6 (1) 

where e is a control parameter. The value of e c = 
corresponds to the saddle-node (tangential) bifurcation 
when the stable and unstable points touch each other in 
x = and disappear. 

Below the critical parameter value (i.e., for e < e c ), 
the stable fixed point is observed, while above e c a nar- 
row corridor between the function f(x) and the bisector 
x n+ i = x n exists, such that the point representing the 
state of the map (pQ) moves along it (Fig. [T]). This move- 
ment corresponds to the laminar phase, its mean length T 
being inversely proportional to the square root of (e — e c ), 
i.e. 

T~(e-e c )- 1 / 2 . (2) 

To develop the theory of type-I intermittency in the 
presence of noise, we consider the same quadratic map (pp) 
with the addition of a stochastic term £ n 

^n+l = %n H~ % n H~ C H~ (3) 

where £ n is supposed to be a delta-correlated white noise 
[(£»> =0, (tntm) =D5(n-m)]. 

The influence of the stochastic term £ n on the behavior 
of the system is governed by the value of parameter D. 
For positive values of the control parameter e (e > 0), the 
point corresponding to the behavior of system ([3]) moves 
in the iteration diagram along the narrow corridor, its 
motion being perturbed by the stochastic force. As far as 
the intensity of the noise is not large, the characteristics 
of classical type-I intermittency are observed. 

A different scenario occurs for control parameters e as- 
suming negative values (e = —e, where e = |e| > 0). 
In this case, the point corresponding to the behavior 
of system (J3j) is localized for a long time in the region 
x < x c = e 1 ! 2 and its dynamics is also perturbed by the 
stochastic force. As soon as the system state point arrives 
at the boundary x c = e 1 ^ 2 due to the influence of noise, 
a turbulent phase arises, though such kind of events is 
very rare. 

In this case, the behavior of the map (pJ]) differs radi- 
cally from the dynamics of the system (pQ), since the tur- 
bulent phases are not observed for e < if there is no 
noise. Therefore, such a region of negative values of the 



e-parameter is the main subject of interest for the type-I 
intermittency in the presence of noise. 

Having supposed that: (i) the value of e is negative 
and rather small and (ii) the value of x changes per one 
iteration insufficiently, we can consider {x n +\ — x n ) as 
the time derivative x and undergo from the system with 
discrete time (|3]) to the flow system, in the same way as in 
the case of the classical theory of the type-I intermittency. 

Since the stochastic term is present in (j3j) we have to 
examine the stochastic differential equation 

dX = (X 2 -e)dt + dW (4) 

(where X(t) is a stochastic process, W(t) a one- 
dimensional Winner process, e = |e|) instead of the or- 
dinary differential equation dx/dt = x 2 + e considered in 
the classical theory of type I intermittency. 

The stochastic differential equation (pfj) is equivalent to 
the Fokker-Plank equation 

dp x (x,t) d 2 Dd 2 p x (x,t) 
^^ = ~^ X - £)pxM) + J^^ (5) 

for the probability density px(x,i) of the stochastic pro- 
cess X(t). Contrarily to what done in Ref. pjl (where 
the backward Fokker-Plank equation was used), we here 
consider the forward one that allows us to obtain the 
explicit form of the distribution of the laminar phase 
lengths. The chosen initial condition is px(x,0) = S(x), 
where S(-) is a delta-function. Such choice of the initial 
form of the probability density px(x,0) corresponds to 
the beginning of the laminar phase, when the point rep- 
resenting the state of the system ([3]) is in the place with 
coordinate x = at time t = 0. In other words, we 
suppose that the reinjection probability is a (^-function 

P in (x)=5(x) (6) 




FIG. 1: (Color online) The iteration diagram for map JTJ) (a) 
e > and (6) e < 0. The stable and unstable fixed points 
of (pp) are shown by • and o, respectively 
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and after the relaminarization process the system is al- 
ways returned to the state x = 0. Although the rein- 
jection probability Pi n {%) is well-known to be important 
factor and should be taken into account when the sta- 
tistical pro perties of the intermittent behavior are stud- 
ied [22I. 123^. in the considered problem the form of the 
reinjection probability practically does not influence on 
the distribution of the laminar phase lengths (and the 
dependence of the mean laminar phase length on the 
criticality parameter, respectively), as it will be shown 
below. 

To reduce the number of the control parameters the 
normalization z = x/y/e, r = ty/e may be used, after 
which Eq. ([5|) may be rewritten in the form 



dp z (z,r) 
dr 



--(( 2 2 -iW(*,r)) + -^A_^ (7) 



where D* = De~ z l 2 , p z (z,r) = p x (^ 1/2 , re' 1 ! 2 ). 

Since the coordinate of the system state is localized 
for a long time in the region z < z c = 1, we suppose 
that the probability density may be written in the form 
Pz(z,t) = A(r)g(z), < 1, where A(r) > decreases 
very slowly as time increases, i.e. dA/dr w 0. The func- 
tion g(z) should also satisfy the conditions 



g(z) > Vz < 1 and 



1 

/ 



](z)dz< 00. (8) 



Under the mentioned assumption, we consider the or- 
dinary differential equation 



D*g"(z)-2((z 2 -i)g(z)y = 

instead of Q for the region z < 1. 
This equation is equivalent to 



(9) 



D*g'(z) - 2(z 2 - l)g(z) + & = 0, (10) 

where C± is constant. To solve this equation we use the 
integrating factor 



2 / 'z 3 

M(z)=exp(-—l--z 



The solution of ([TO]) may be found in the form 

n 3 



(ii) 



Ci 



(12) 



Since the obtained g(z) satisfies the conditions (|BJ) only 
if C\ =0, the probability density pz(z,r) in the region 
z < 1 is 



(13) 




FIG. 2: (Color online) The evolution of the probabil- 
ity density px(x,t) obtained by means of the direct nu- 
merical integration of Fokker-Plank equation d5}, £ — 10 -2 , 
D — 2.5 x 10 -4 . (a) The initial fragment of the density evo- 
lution involving the transient (0 < t < tt r , Ur ~ 30). (b) The 
long-time evolution of px(x, £), with the transient being omit- 
ted, t > 50 



The decrease of A(r) should be determined by the 
probability distribution taken in the boundary point 
z = 1, i.e., dA(r)/dr ~ —pz(l,r). This assumption 
may be rewritten as 



- -kA(r)ex.p ^ 



dr 



3L>* 



(14) 



where k is a proportionality coefficient. Evidently, the 
decrease of A(r) is described by the exponential law 

A(t) = A(0) exp(-^r), rj = exp (-4/(3L>*)). (15) 

Having returned to the initial variables x and t we 
obtain the following expression for the probability density 
Px(x,t) 



where 



and 



A(t) = A(0) exp 



1 /4e 3 / 2 \ 



(16) 



(17) 



(18) 



with A(t) being considered as a normalizing factor, i.e., 



( ex - ^- ) 



To confirm the assumptions made above and the ob- 
tained equations, we have compared the evolution of the 
probability density px(x,i) given by (fl6|) with the re- 
sult of the direct numerical calculation of the Fokker- 
Plank equation ([5]) with the values of control parameters 
e = 10~ 2 , L> = 2.5 x 10~ 4 . 

The evolution of the probability density px(x,t) 
obtined by the numerical calculation of ([5]) is shown in 
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FIG. 3: (Color online) The profiles of the probability density 
px(x,t*) taken in the different moments of time t* obtained 
from the direct numerical calculation of Fokker-Plank equa- 
tion ©. 

Fig. [21 One can see that after the very short transient 
< t < t tr the probability density px(x^t) arrives the 
state being close to stationary (Fig. [2]a). After that 
the value of px(x^t) decreases very slowly (according 
to the exponential law) with time increasing, with the 
form of the dependence of the probability density on x- 
coordinate being invariable (Fig. [2]&). 

Fig. [3] also shows the profiles of the probability den- 
sity px{%-,~k*) taken in the different moments of time. 
It is evident, that after a very short transient (curve 1, 
t* = 10), the density px(x^t) practically does not change 
when time increases. 

Two different profiles px(x,t*) corresponding to the 
time moments t% = 3 x 10 1 and ^ = 2x 10 4 (curves 2 
and 3, respectively) are very close to each other despite of 
the very large time interval At = t\ — t\ between them. 
Moreover, they are in very good agreement with the ap- 
proximated solution A(0)g(x) described by Eq. ([T6]) and 
shown in Fig.[3]by means of squares. As time goes on, the 
amplitude of the probability density decreases accord- 
ing to the exponential law, but very slowly (see Fig. O 
curves 4 and 5, t% = 2.5 x 10 5 and t% = 10 6 , respectively), 
although the probability density form remains the same 
for all times. 

Therefore, taking into account the results of the di- 
rect numerical calculations of Fokker-Plank equation ([5]) 
and the comparison with the obtained approximated so- 
lution ([T6]h we come to the conclusion that our assump- 
tions are correct and can be used for the further analysis. 

The evolution of the probability density px(%,t) may 
be considered separately on two time intervals < t < t tr 
and t tr < t < +oo, respectively. The first time interval 
corresponds to the transient when the probability den- 
sity px(x,t) evolves to the form ([T6|) being close to sta- 
tionary. Only when < t < t tr the form of the reinjec- 
tion probability Pi n (x) may influence on the evolution 
of the probability density px(x,t). For t > t tr (when 
the transient is elapsed), the evolution of the probabil- 
ity density is defined completely by Eq. ([T6]) and it does 



not depend entirely on the reinjection probability Pi n (x). 
Since the transient is very short in comparison with the 
exponential decrease of the probability density px(%,t) 
we can neglect them and use only the second time interval 
t tr <t < +oo to obtain the statistical characteristics of 
the type-I intermittent behavior in the presence of noise. 
It is clear, that in this case the obtained results do not de- 
pend on the relaminarization process and the reinjection 
probability P in (x) . 

The distribution p(t) of the laminar phase lengths t 
may be defined from the relationship between px(%,t) 
and p(t) 

p(t) = -J 9 -^d X . (20) 

-co 

Using relations ([T6jh (|T7|) and ([T9|) one can obtain, that 
the laminar phase distribution is governed by the expo- 
nential law 

p(t) = T- 1 exp(-t/T), (21) 

where T defined by Eq. ([T8]) is the mean length of the 
laminar phases. The obtained expression ([T8]) for the 
mean length T of the laminar phases is consistent with 
the formal solution T ~ e~ 1 ^ 2 f(a 2 s~ 3 ^ 2 ) derived in the 
previous considerations [I?], [y| [24| and coincides with 
the formula of T given in [19[. If the criticality pa- 
rameter e is large enough, the approximate equation 
InT ~ Z) _1 £ 3 / 2 may be used (see [19] for detail). Based 
on the consideration carried out above we state that 
Eqs. (jTHJ) and ([2T]) do not depend practically on the re- 
laminarization process properties and may be used for 
the arbitrary reinjection probability Pi n {x). 

II. SAMPLE SYSTEM DYNAMICS 

To verify the obtained theoretical predictions, we con- 
sider numerically two different dynamical systems show- 
ing type-I intermittency, with a stochastic force being 
added. As such test systems we have selected (i) the 
quadratic map and (ii) driven Van der Pol oscillator. 

A. Quadratic map with stochastic force 

We start considering the interplay between type-I in- 
termittency and noise using the quadratic map 

+ A + e + L>£ n , modi, (22) 

where the "mod 1" operation is used to provide the re- 
turn of the system in the vicinity of the point x = 0; 
A = 0.25, and the probability density of the stochas- 
tic variable £ is distributed uniformly on the interval 
£ G [—1, 1]. The map ([22]) may be brought to (J3j) with the 
help of a linear variable transformation. If the intensity 
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FIG. 4: (Color online) (a) The dependence of the mean length 
T of the laminar phases on the criticality parameter e = \e\ 
(e < 0) for ([22]). The points obtained by the iteration of (12^1) 
are shown by symbols (•). The theoretical law InT ~ e 3 ^ 2 is 
shown by the solid line, (b) The distribution of the laminar 
phase lengths for map ([22]) . the criticality parameter value 
has been selected as e = 10 -11 . The theoretical exponential 
law ([21[) is shown by the solid line 



of noise D is equal to zero the saddle-node bifurcation 
is observed for e = 0. The type-I intermittent behav- 
ior is observed for e > 0, whereas the stable fixed point 
takes place for e < 0. Having added the stochastic force 
D = 10 -7 in ([22]) we can expect that the intermittent be- 
havior may be also observed in the area of the negative 
values of the criticality parameter e. 

The dependence of the mean laminar phase length T 
on the criticality parameter below the point e c = is 
shown in Fig. [4j To compare it with the obtained theo- 
retical prediction InT ~ D~ 1 e 3 ^ 2 the abscissa in Fig.0] a 
has been selected in the £ 3 / 2 -scale [e = |e|, e < 0), 
whereas the ordinate axis is shown in the logarithmic 
scale. One can see the excellent agreement between the- 
oretical law ([T8]) and data of numerical calculation. 

The distribution of the laminar phase lengths is also 
in very good accordance with the exponential law ([2T]) 
predicted by the theory of the type-I intermittency with 
noise (see Fig. [4] b). Note the presence of the small re- 
gion of the short laminar phase lengths in Fig. 2J b where 
the deviation from the prescribed exponential law ([2T]) 
is observed. This region corresponds to the transient 
< t < t tr when the probability density px(x,i) evolves 
to the form ([T6]) being close to stationary as it was dis- 
cussed in Sec. |T] The existence of this transient time 
interval does not influence practically on the character- 
istics ([T8]) and ([2T]) of the intermittent behavior in the 
presence of noise in the full agreement with the conclu- 
sions made above. 

So, the intermittent behavior observed in the quadratic 
map with the stochastic force agrees well with the theo- 
retical predictions obtained in Sec. H Since the theory of 
type-I intermittency has been developed on the basis of 
the model ([3]) being very close to map ([22]) , it is manda- 
tory to examine another system to ensure that our theo- 
retical conclusions are correct and applicable for a wide 



spectrum of nonlinear systems. 

B. Van der Pol oscillator driven by the external 
harmonic signal in the presence of noise 

We consider as a second model the system given by a 
van der Pol oscillator 

x - (A - x 2 )x + x = A sin(uj e t) + D£(t) (23) 

driven by an external harmonic signal with the ampli- 
tude A and frequency uo ei with an added stochastic term 
D£(t). The values of the control parameters are selected 
to be A = 0.1, oo e = 0.98. For these control parameters 
and for D = 0, the dynamics of the driven van der Pol 
oscillator becomes synchronized when A = A c = 0.0238. 
The probability density of the random variable is 

m = vk^(-&)- (24) 

where a 2 = 1. To integrate Eq. ([23]) the one-step Euler 
method has been used with time step h = 5 x 10 -4 , the 
value of the noise intensity has been fixed as D = 1. 

It is well-known that under certain conditions (i.e., for 
the periodically forced weakly nonlinear isochronous os- 
cillator), the complex amplitude method may be used to 
find the solution describing the behavior of oscillator ([23]) 
without noise in the form x(t) = Rea(t)e zu;t . For the 
complex amplitude a(t) one obtains the averaged (trun- 
cated) equation a = — iva + a — \a\ 2 a — ik, where v is the 
frequency mismatch, and k is the (re- normalized) ampli- 
tude of the external force. For the small v and large k, 
the stable fixed point on the plane of the complex ampli- 
tude a* = const corresponds to the synchronous regime, 
with the synchronization destruction corresponding to 
the saddle-node bifurcation associated with the global 
bifurcation of the limit cycle birth [25|, [26[. Therefore, 
below the boundary of the synchronization regime (for 
small values of the frequency mistuning), the dynamics 
of the phase difference 

Atp(t) = cp(t) - uo e t (25) 

(where cp(t) is the phase of the driven oscillator) demon- 
strates time intervals of phase synchronized motion (lam- 
inar phases) persistently and intermittently interrupted 
by phase slips (turbulent phases) during which the value 
of \Aip(t)\ jumps up by 2tt. The mean length T of the 
laminar (synchronous) phases depends on the criticality 
parameter e = (A c — A) according to the power law ([2]) 
corresponding to the type-I intermittency. 

If the stochastic term D£(t) is added (D ^ 0) the man- 
ifestation of the regularities of type-I intermittency with 
noise described above is revealed in the parameter range 
A > A c (see Fig.[5j a). For the negative values of the crit- 
icality parameter e the law InT ~ D~ 1 e 3 ^ 2 is expected to 
be observed. To make this law evident, the abscissa in 
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FIG. 5: (Color online) (a) The dependencies of the mean 
length T of the laminar phases on the criticality parameter 
e — (A c — A) for the driven van der Pol oscillator with the 
stochastic force (|23p . The points obtained by the numeri- 
cal integration of ()23|) are shown by symbols •. The theo- 
retical law lnT^s 3 / 2 is shown by the solid line, (b) The 
distribution of the laminar phase lengths for the driven van 
der Pol oscillator, with the amplitude of the external signal 
A = 0.0245 being taken above the critical point A c = 0.0238 
(e = -7 X 10" 4 ). The ordinate axis is presented in the loga- 
rithmic scale. The theoretical exponential law (|21[) is shown 
by the solid line 



Fig. [5j a has been selected in the £ 3 / 2 -scale (e = |e|) and 
the ordinate axis T is shown in the logarithmic scale. One 
can see again the excellent agreement between the numer- 
ically calculated data and theoretical prediction ([T5|). 

The distribution of the lengths of the laminar phases 
N(t) obtained for A > A c also confirms the theoretical 
results given above. Indeed, the distribution N(t) shown 
in Fig. O b is in the very good accordance with the the- 
oretically predicted exponential law (f2T]h with the small 
region of the short laminar phase lengths deviating from 
the exponential law being revealed as well as for the dis- 
crete map ([22]) that corresponds to the short transient 
taking place after the relaminarization. Again, as for the 
discrete map considered in Sec. Ill Al the existence of this 
transient time interval does not distort the characteris- 
tics of the intermittent behavior observed in the presence 
of noise. 



III. EXPERIMENTAL OBSERVATION OF THE 
CHARACTERISTICS OF TYPE-I 
INTERMITTENCY IN THE PRESENCE OF 
NOISE 

In parallel with the numerical analysis of the type-I 
intermittent behavior with noise we have also studied 
experimentally the dynamics of the periodic oscillator 
driven by the the external harmonic signal in the pres- 
ence of noise to confirm the theoretical and numerical 
results given in Sec. UlandHIl In the experiment we have 
used the simple electronic oscillator where all parameters 
(including noise amplitude) may be controlled precisely. 

The experimental setup is shown in Fig. [6j The ba- 
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FIG. 6: The schematic diagram of the experimental setup. 
The control parameters have been selected as the following: 
Ri = 10 Ohm, R 2 = 630 Ohm, R 3 = 56 Ohm, L 3.3 niH. 
C = 150 nF, Ci = 330 nF. The operational amplifiers OP1 
and OP2 are both of the TL082 type and the operational 
amplifier OP3 is of the TDA2030 type 



OP1 



OP2 



input 




FIG. 7: The schematic diagram of the nonlinear converter. 
The control parameters have been selected as the following: 
Ri 2.7 kOhm, R 2 7.5 kOhm, R 3 100 Ohm, R 4 
7.5 kOhm, R 5 = 12 kOhm, R 6 = 4.7 kOhm. The diodes Di 
and D 2 are of the 1N4148 type. The operational amplifiers 
OP1 and OP2 are both of the TL082 type and the operational 
amplifier OP3 is of the LF356N type 



sis element of the scheme we use the generator with the 
linear feedback and nonlinear converter (NC) [27]. The 
diagram of the nonlinear converter is shown in Fig. 
The characteristics of nonlinear converter were controlled 
with resistor Rq (see Fig. [7J). Since the generator demon- 
strates both chaotic and periodic oscillations, the control 
parameters of it have been selected in such a way for the 
generated signal to be periodic. The frequency of the 
autonomous periodic oscillations was 8.805 kHz. As a 
source of driving harmonic signal the MOTECH-FG503 
functional generator (FG) has been used. The behav- 
ior of the oscillator driven by the external harmonic sig- 
nal in the presence of noise was analyzed by means of 
the Agilent E4402B spectrum analyzer and L-Card L- 
783 analog-digital converter (ADC) PCI-card with 12-bit 
resolution. 

The noise generator (NG) shown in Fig. [HI a provides 
the noise signal being close to Gaussian one [28|. The 
distribution of noise p(V) is shown in Fig. [HI b. The in- 
tensity of noise may be controlled easily by means of the 
variation of the potentiometer R4. 

As well as for Van der Pol oscillator driven by the exter- 
nal harmonic signal in the presence of noise (see Sec. Ill B| ) 
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FIG. 8: (Color online) (a) The schematic diagram of the noise 
generator. The control parameters have been selected as the 
following: Ri = 72 Ohm, R 2 = 1 kOhm, R 3 = 1 kOhm, 
R A 10 kOhm, R 5 10 kOhm, Ci 1 //F, C 2 0.1 /xF, 
C 3 = 470 /xF, *7 = 15 V. The diode D x is of the BZX79-C11 
type. The operational amplifier OP1 is of the TL082 type, 
the transistor T\ is of the BC546B type. (6) The distribution 
p{V) of noise 



10 4 



10 



10 2 



N(t) 



,3- 



.* 



io 1 - 



io u 



6.0x10" 



1.2x10- 



t, [s] 

1 

1.8x10- 



below the boundary of the synchronization regime (for 
small values of the frequency mistuning) the dynamics of 
the difference 

A<p(t) = <p(t) - <p e (t) (26) 

between the phase of the driven oscillator cp{t) and 
the phase of the external harmonic signal ip e (t) should 
demonstrate time intervals of phase synchronized mo- 
tion (laminar phases) persistently and intermittently in- 
terrupted by phase slips (turbulent phases) during which 
the value of \ Acp(t)\ jumps up by 2tt. The distribution of 
the laminar phase length is expected to obey to the ex- 
ponential law ([2T]) predicted by the theory of the type-I 
intermittency with noise. 

Since the dependence of the mean laminar phase length 
on the criticality parameter has already been studied ex- 
perimentally 20} in our experiment we focus on the con- 
sideration of distribution of the laminar phase lengths. 
These distributions N(t) obtained experimentally for the 
different values of the amplitude V m of the external har- 
monic signal and noise intensity D are shown in Fig. [9j 
The frequency of the external harmonic signal has been 
fixed as / = 8.75 kHz, the values of the amplitude V of 
the external signal have been selected in such a way for 
the driven oscillator to be synchronized if the intensity of 
noise is equal to zero, i.e., Vm > V c , where V c = 145 mV 
is the amplitude of the external signal corresponding to 
the synchronization threshold. In the presence of noise 
the phase slips are revealed and the intermittent behav- 
ior is observed. One can see that the distributions of the 
lengths of the laminar phases N(t) are in the very good 
accordance with the theoretically predicted exponential 
law ([2T]) . The small regions of the short laminar phase 
lengths where the deviation from the prescribed exponen- 
tial law ([2T]) is observed also take place as well as in the 
case of the numerical simulations of the model systems 
considered in Sec. [Ill Therefore, we come to conclusion 
that the experimental observations confirm the obtained 
theoretical results concerning the type-I intermittent be- 
havior in the presence of noise. 



FIG. 9: (Color online) The distributions of the laminar 
phase lengths for the driven periodic oscillator in the pres- 
ence of noise obtained experimentally. The amplitude V m 
of the external signal and the noise dispersion a have been 
selected the following: (♦) V m = 170 mV, a = 475.17 mV; 
(■) V m 170 mV, a 298.11 mV; (•) V m 150 mV, 
a = 141.72 mV. The ordinate axis is presented in the loga- 
rithmic scale. The approximations corresponding to the the- 
oretical exponential law (|21[) are shown by the dashed lines 



IV. CONCLUSIONS 

In conclusion, we have reported a type of intermit- 
tency behavior caused by the cooperation between the 
deterministic mechanisms and random dynamics. Hav- 
ing examined three sample systems both numerically and 
experimentally we can conclude that (i) noise induces 
new features in the intermittent behavior of a system 
demonstrating type-I intermittency, with new dynamical 
properties being observed above the former value of the 
criticality parameter; (ii) the results of numerical sim- 
ulations and experimental observations are in excellent 
agreement with the developed theory; (iii) the statistical 
characteristics of the perturbations of type-I intermit- 
tency as well as the relaminarization process properties 
and the reinject ion probability do not seem to play a ma- 
jor role. Though the characterization of the intermittent 
process has been explicitly derived here for model sys- 
tems, we expect that the very same mechanism can be 
observed in many other relevant circumstances where the 
level of natural noise is sufficient, e.g. in the physiologi- 
cal [29|, [30[ or physical systems [H . 
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